Gravitational collapse in non-minimally coupled gravity: 
finite density singularities and the breaking of the no-hair theorem 
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In this work we study the dynamics of gravitational collapse of a homogeneous dust sphere in a 
model exhibiting a linear non-minimal coupling between matter and curvature. The evolution of 
the scale factor and the matter density is obtained for different choices of Lagrangean density of 
matter, highlighting the direct physical relevance of the latter in this theory. Following a discussion 
of the junction conditions and boundary terms in the action functional, the matching with the outer 
metric and event horizon are analyzed. 

We find that a distinct phenomenology arises when compared with standard results for the 
Oppenheimer-Snyder collapse, namely the possibility of finite density black holes and the break- 
ing of the no-hair theorem, due to a dependence of the end state of a black hole on the initial radius 
of the spherical body. 

PACS numbers: 97.60.Lf,04.20.Fy, 98.35.Cc 



I. INTRODUCTION 

Two of the major challenges faced by contemporary 
cosmology arc the nature of dark energy [1, 2] and dark 
matter [3] (or perhaps a unification of the two [4], which 
account for ~ 96% of the Universe, and to gather some 
insight into a more encompassing theory of gravity. A 
rather straightforward approach for the latter problem 
resorts to the substitution of the linear scalar curva- 
ture term in the Einstein-Hilbert action with a function 
of the scalar curvature, f(R) [5] or other scalar invari- 
ants of the theory: extensions relying on a functional 
dependence of the action on the Gauss-Bonnet invari- 
ant G = R 2 - AR ltu R> iV + R a p^R a0fnJ are the most well 
studied models, given their invoked origin in a low-energy 
effective description of String Theory and strong impli- 
cations in braneworld scenarios [6]. 

The more tractable class of f(R) theories have had con- 
siderable success in replicating the early period of rapid 
expansion of the universe, as shown by the Starobinsky 
inflationary model f(R) = R + aR 2 [7]. At late times, 
the accelerated expansion of the Universe has also been 
addressed suitably [8]. Solar system tests, mostly arising 
from the parameterized post-Newtonian (PPN) metric 
coefficients derived from this extension of General Rela- 
tivity (GR), have also been discussed [9]. A clear phe- 
nomenological consequence of f(R) theories it the addi- 
tion of an increasing, repulsive contribution to the New- 
tonian potential, for power law terms f(R) = foR n [10]. 
Aside from the more usual metric affine connection (that 
is, where the affine connection is taken a priori as de- 
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pending on the metric), the so-called Palatini approach 
[11] (where both the metric and the affine connection are 
taken as independent variables) has also been considered. 

Further expanding on this elegant generalization of 
GR, another interesting possibility has arisen: that the 
coupling between matter and geometry in the Einstein- 
Hilbert action is non- minimal — i.e. not enforced solely 
by the invariant y rr gd 4 x and the use of the metric to 
raise and lower indexes and the associated covariant 
derivative. A non-minimal coupling would imply that 
geometric quantities (such as the scalar invariants) could 
explicitly appear in the action [12] (see also Ref. [13] 
for an early proposal in the context of Riemann-Cartan 
geometry) . 

Indeed, one is motivated to do so by the presence of a 
non-minimal coupling stemming from one-loop vacuum- 
polarization effects in the formulation of Quantum Elec- 
trodynamics in a curved spacetime [14], as well as in the 
context of scalar-tensor theories, when considering mat- 
ter scalar fields [15]. Furthermore, it has been shown that 
a non-minimally coupled model cannot follow the usual 
procedure establishing the equivalence between f(R) and 
a single scalar-tensor theory [16]; indeed, while a model 
with two scalar fields may describe the same dynamics as 
Eq. (1), it still requires one of these fields to appear non- 
minimally coupled with the matter Lagrangean density 
[16-20]. 

The presence of non-minimal coupling leads to sev- 
eral phcnomcnological consequences: in particular, it im- 
plies that regions with extreme curvature could lead to 
considerable deviations from the dynamics predicted by 
Einstein's theory [12]. A wide range of results has un- 
folded, including the impact on solar observables [21], 
axisymmetric astrophysical solutions [22] , the possibility 
to account for galactic [23] and cluster dark matter [24] , 
a mechanism for mimicking a Cosmological Constant 
at astrophysical scales [17], post- inflationary reheating 
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[25] or the current accelerated expansion of the universe 
[26] (also including the so-called "extended quintessence" 
[27]). Finally, a thorough discussion of the relevance of 
the choice of Lagrangean density for a perfect fluid and 
its direct impact on physical observables was discussed 
in Ref. [28]. This choice will be of the utmost relevance 
in the current work. 

In this study, one addresses how a non-minimal cou- 
pling modifies gravitational collapse; similar studies have 
been performed in the case of standard f(R) theories (i.e. 
with = f(R) ^ R and f 2 (R) = 1) [29-31]. Sev- 

eral simplifications are made, namely that the collapsing 
body is purely spherical body composed of a homoge- 
neous distribution of dust — similarly to the well-known 
Oppenheimer-Snyder (OS) collapse [32] thoroughly stud- 
ied in GR. Also, the simplest, linear non-minimal cou- 
pling f2{R) ~ R is considered, and a trivial curvature 
term fi (R) — R is taken so to highlight the effect of the 
former on the collapse. 



II. THE MODEL 

One considers a model that exhibits a non-minimal 
coupling between geometry and matter, as expressed in 
the action functional [12], 



[nf 1 (R) + f 2 (R)C]V^d 4 x 



(1) 



Variation with respect to the action yields the modified 
Einstein field equations, 



1 

n 



F\ + —F 2 C J G^u — —jiTfiv 
Zk 



(2) 



A M „ I F 1 + -F 2 c] + \ 9llv (h - FiR - -F 2 RC 



with k = c /16irG and = V M V„— g M i,Q The energy- 
momentum tensor is defined as 



J2 



(5) 



Again, in the absence of a non-minimal coupling, 
f2(R) = 1 and the covariant conservation of the energy- 
momentum tensor is recovered. This feature implies that 
the motion of the matter distribution described by a 
Lagrangian density C does not follow a geodesic curve. 
Clearly, a violation of the Equivalence Principle may 
emerge if the r.h.s. of the last equation varies signifi- 
cantly for different matter distributions, which suggests 
a method of testing the model and imposing constraints 
on the associated couplings. This feature is a fundamen- 
tal characteristic of a non-minimally coupled model, as 
shown in Ref. [18]. 



III. GRAVITATIONAL COLLAPSE OF A 
HOMOGENEOUS FLUID 

A. Linear non-minimal coupling 

In Ref. [28] , it was argued that the correct Lagrangean 
density C for a perfect fluid is C = —p, as the non- 
minimal coupling disables the usual on-shell equivalence 
with other forms (such as C — p). Notwithstanding, in 
this study one aims at addressing both forms, so that 
the impact of choosing a particular Lagrangean density 
on gravitational collapse may be gauged directly. 

Since one considers a homogeneous dust distribution, 
i.e. p = const., p — 0, this is attained by writing C = 
—ap, with a = — > C = p and a = 1 — > C = —p. A dust 
distribution implies that there is no supporting pressure 
to prevent collapse, with no shell crossing during the later 
(and no exchange of momentum). 

Inserting the above expression, together with a triv- 
ial curvature term fi(R) = R and a linear non-minimal 
coupling 



f 2 (R) = l + -R 



(0) 



-9 Sg^ 



(3) and the energy-momentum tensor for a dust distribution, 



As expected, GR is recovered from Eq. (3) by setting 
f x (R) = R and f 2 (R) = 1. 
The trace of Eq. (3) reads 



F\ + -F 2 C) R = i/ 2 T - 3D (f 1 + -F 2 C) + 2f 1 , 

k / Zk \ k J 

(4) 

where T is the trace of the energy-momentum tensor. 

Resorting to the Bianchi identities, one concludes that 
the energy-momentum tensor of matter may not be co- 
variant ly conserved, as 



T^v = pu^ 



(7) 



where the velocity obeys u^u v = — 1 and u M V ' = 0, 
one finds that Eq. (3) reads 



( 



1 ea \ n 
1 2 P } <J > 1/ 



(8) 



1 / e \ lea eot . 

Its trace yields 
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V. EVOLUTION OF THE GRAVITATIONAL 
COLLAPSE 

= k P - 3eoQp 
2 K 2 + e(2a-l)p ' 

One now derives the equations of motion driving the 
dynamics of gravitational collapse. With the assumed 
The non-conservation of the energy-momentum tensor FRW metric, the scalar curvature given by Eq. (9) bc- 
Eq. (5) becomes comes 



eR 



(10) 



IV. ADAPTABILITY OF THE FRW METRIC 

In GR, the study of OS collapse is tantamount to the 
determination of the evolution of the scale factor a(t) ap- 
pearing in the Friedmann- Robertson- Walker (FRW) met- 
ric, as given by the line element 



R = 



np + 3ea (p + 3§p) 
2k 2 + e(2a-l)p 



(14) 



Thus, the temporal component of Eq. (10) becomes 



Pit) 



Po 



a 

3- + — 

a k 

K + eR 



-eR 

a-l 



(1 -a)R 



(15) 



f-V 

eRo J V a J 



ds 2 



-dt 2 +a 2 (t) 



dr 2 



1 - kr 2 



+ r 2 {d6 2 + sin 2 dd(f> 2 ) 



(11) 

where k is the (negative) spatial curvature. 

The use of the above metric implies that one has a 
position-independent scalar curvature R: naturally, this 
stems from the identification R — p(t) that arises from 
the trace of the unperturbed Einstein equations. Naively, 
one could expect that in the current scenario a homoge- 
neous density also gives rise to a scalar curvature ex- 
hibiting only a time dependence, albeit a more convo- 
luted one, as seen in Eq. (9). This, however, fails to 
acknowledge that a more general metric could give rise 
to a radial dependence of R, via the terms involving the 
(space dependent) metric components appearing in the 
D'Alembertian operator Q 

To assess this, one therefore adopts a more general 
metric, given by 



ds 2 = -dt 2 + U(r,t)dr 2 
U{r,t) = A\{t)h{r) , 



-V(r,t)(d6 2 +sin 2 
V(r,t) = A 2 2 (t)r 2 



(12) 



so that t measures time according to a in-falling observer 
(see Ref. [29] for a similar treatment in f{R) theories). 
One computes the required differential operator, 



(13) 



which, by construction, is a function of time only. Thus, 
from Eq. (9), it becomes clear that the scalar curvature 
indeed inherits the homogeneity of the matter density 
distribution, R = R(t), and the FRW metric is indeed 
admissible in the present scenario. 



having used = (1,0) in the adopted comoving coor- 
dinates, and where o denotes initial values (without loss 
of generality, one henceforth sets ao = a (to) = 1). In- 
spection shows that the radial component of Eq. (10) is 
trivially null. 

Substituting Eq. (14) into the above and considering 
that a — 0, 1 is a "binary" variable yields, after some 
algebra, the closed expression 



P(t) 



p a 



l + e(l 



a) -22 



2h- 



(a-3-1) 



(16) 



Even without the knowledge of the evolution of the 
scale factor a(t), this result is revealing of the effect of 
the adopted non-minimal coupling: if C = —p — > a = 
1, it leads to an unchanged evolution for the increasing 
density, p ~ a~ 3 . 

If, however, C = p — s- a = 0, one finds that as the 
spherical body collapses and a(t) — > 0, it is not infinitely 
compressed, p — > oo, but rather it tends towards a fi- 
nal state of finite density Pf = 2k 2 /e ~ Mp/e, where 
Alp is the Planck mass. Furthermore, if e ^> 1, this 
density is low enough that no (yet unknown) quantum 
effects should arise and one is left with a classical object 
of vanishing size but finite density! 

Such counter-intuitive result is of course related with 
the non-conservation of energy, the most striking feature 
of non-minimally coupled models. Also, notice that the 
choice of Lagrangean density indeed has a crucial role in 
the physical outcome of the model, as discussed in the 
previous section. 

Replacing Eq. (14) into the modified field Eqs. (9) and 
following some tedious computations, one obtains the two 
dynamical equations of motion relating the density with 
the scale factor and the spatial curvature, 



= (k 2 + 2eap) (k + 2ad) + (k 2 - Aeap) {of 
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-na 2 p = (k 2 ~ep)k + [k 2 + (3a - l)ep] (a) 2 + 
e(a — Vja'dp . 



(17) 



Notice that a = leads to the presence of the sec- 
ond derivative of the scale factor in the second equa- 
tion. Evaluating the first one at t — to, and assuming 
a collapse with initial null velocity, <z(0) = 0, one gets 
k = — 2d > 0, regardless of a (with a = a(0)) — unless 
a = 1 and po = —K 2 /2e (implying a negative coupling 
strength e), in which case k remains undefined. 
One may use Eq. (17) to eliminate d and write 



one finds that Eq. (22) becomes the usual evolution equa- 
tion found in GR, 



a'(rj) = -yfa{\- a) , o(0) = 1. 



(25) 



The effect of the non-minimal coupling is ascribed to the 
relation between the cycloid time r\ and the coordinate 
time t: if C = p — > a = 0, this amounts only to a shift 
of the spatial curvature k; if C = —p — > a = 1, Eq. (24) 
leads to a more complex relation between the latter. 
As in GR, the solution of Eq. (25) is 



-na 2 p = [2k 2 - (1 + a)ep] k + [2k 2 + (5a - l)ep] (a) 
which, again assuming a(0) = 0, yields 



2 

(18) 



1 + cos(ry) 



(26) 



so that the gravitational collapse ends when r/ = tt. 



k = 



Pi i 



3K2-(l + a)e ' 



(19) 



defining the dimensionless parameter eq = epa/n 2 . Thus, 
one finds that the non-minimal coupling induces a shift 
from the value for the spatial curvature found in GR, 
kg = pq/6k. Moreover, a positive spatial curvature also 
implies that the initial density and the coupling strength 
are constrained by 



Po < 



e (l + o) 



(20) 



Substituting the expression for pit) found in Eq. (16) 
into the first of Eqs. (17) and using the identity a(a — 
1) = 0, one obtains 

= (k 2 + 2eap a~ 3 ) (k + 2aa) + (k 2 - 4eap a~ 3 ) (a) 2 , 

(21) 

leading to a simplified equation of motion for the scale 
factor 



, n a 2 + (a + l)ae 
a — — \ I k(l — a) ^ 



cr + 2aeo 



(22) 



supplemented by the initial condition a(0) = 1. 

One may easily check that when a ~ the scalar cur- 
vature behaves as 



R ~ (1 + 2a) 



3fc 



(23) 



so that both choices of Lagrangcan density lead to a sin- 
gularity. 

Defining the cycloid parameter rj as, 



Vk / a 3 + a(a + l)ae 

ai] = — \ — at 

a W a J + zaeo 



(24) 



VI. BOUNDARY MATCHING 

Assuming that space outside the collapsing spherical 
body is empty, Eq. (3) become the trivial vacuum equa- 
tions Rfn, = found in GR, and as such spacetime is 
described by a Schwarzschild metric. The outer region is 
endowed with a coordinate system that does not coincide 
with the one used in the interior region: from spherical 
symmetry, one sees that only the time and radial coor- 
dinates will be different. Labeling these as t 1 and r', 
respectively, one writes the outer metric via the usual 
line element 



ds 2 = - 1 - 



2GM\ ,, 2 
> dt' 2 



2GM\ 
1-— J 



-l 



(27) 



dr 
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r' 2 (d9 2 



For the full spacetime to be well defined, one must 
ensure that the FRW metric valid Eq. (11) inside 
the collapsing body matches smoothly with this outer 
Schwarzschild metric. 

One first requires that the induced metric on the 
boundary is equal on both sides, as any discontinuities 
would lead to an ill-defined scalar curvature R — thus 
formulating the so-called first junction condition. For 
this, one first defines the boundary of the collapsing 
spherical body as a spacelike hypcrsurface given by the 
condition r = 7'* = const, (since r is a comoving coordi- 
nate) or r' = R*(t') (i.e. an external observer sees the 
boundary receding towards R* =0). 

The tangent vectors are given by 



dy a 
dx a 



(28) 



where y a are four coordinates used in the inner and outer 
metrics and x a (a = t, 9, (f) are the three coordinates 
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parameterizing the boundary [33]. The induced metric 
hat, on this hypersurface, defined as, 

Kb = g a ^ a a 4> , (29) 
is then given by the line element 



1 = 0. In the interior of the spherical body I < and 
H~(l) = 1, while outside H~{1) = (and conversely for 
H+). 

The computation of the Ricci tensor and the scalar 
curvature involve derivatives of the metric; recalling that 
H'(l) = 8(1), one obtains 



ds^ = -dr 



a 2 rl(d9 2 
2GM 



Rl{d0 2 + sin 2 



2GM\ 



R* J 



dr' 

IF 



(30) 
dt' 2 + 



R = H~R- +H + R* 



SA 



(33) 



By inspection, one obtains the matching conditions for 
the inner and outer coordinate systems, 



R* 



(31) 



V = 



R*[R* — 2GM + R*(R*) 2 ] 



R* - 2GM 



where the dot still indicates differentiation with respect 
to t. 

The smooth crossover between the inner and outer de- 
scription of spacetime also demands that the derivatives 
of the corresponding metrics are properly matched. A 
quick and dirty approach is to require that the derivatives 
with respect to the radial coordinates are equal; how- 
ever, this procedure is explicitly coordinate-dependent. 
A more elegant, coordinate invariant approach relies on 
the computation of the (dis) continuity conditions for the 
extrinsic curvature tensor K ao [33]. 

In GR, this ensuing second junction condition (in the 
absence of a thin shell, see Ref. [34] for a related dis- 
cussion) translates into the equality K~ b = K~^ b . Erro- 
neously, some assume that this is a universal statement 
born from geometrical necessity or aesthetic considera- 
tions: for example, in Ref. [30], the unwarranted use of 
this equality in the context of f(R) theories leads to a 
time dependent mass M(t'), thus producing an unphysi- 
cal outer vacuum with a non- vanishing, time dependent 
curvature! 

Given the above, one concludes that when GR is mod- 
ified, the second junction condition may in general read 
[K a b] = K^b ~ Kab ^ [35]; with the r.h.s. reflecting the 
altered structure of the equations of motion. 

In order to pursue its correct formulation, two equiv- 
alent paths are available, as described schematically be- 
low. The first works at the level of the equations of mo- 
tion and resorts to a functional description of the relevant 
quantities [33], thus writing the metric as 



g flv =H-(l)g-+H + (l)g+ 



(32) 



where -ff ± is the Heaviside step function and / is an affinc 
parameter describing the crossing of the hypersurface at 



where A = A(g^ v ,g^ v ^ a ,np) and np is the unit normal 
to the hypersurface; the latter enables the relation be- 
tween the induced and the inner and outer metrics at 
the boundary, 

g°P = n a vP + h ah e a a el . (34) 

Following the same procedure for the matter content, 
one writes its energy-momentum tensor as 



T^ = H-(l)T-+H + (l)T+ 



(35) 



with an additional 8(1) term if a boundary layer is 
present. 

By inserting the above expansions for the relevant 
quantities into Eq. (3) and demanding that terms in 
6(1) vanish, one is in principle able to obtain the desired 
second junction condition. Since many modifications of 
GR increase the order of the differential operators in the 
equations of motion, in general this will lead to added 
terms in 6 that yield a discontinuity of the extrinsic cur- 
vature across the hypersurface. Also, particular care 
should be taken with the appearance of crossed terms 
such as H~6, as these are ill-defined as functionals - 
an issue that could be surpassed with the substitution of 
the Heaviside step and Dirac delta functions by suitable, 
convergent, approximations. 

The second, much more elegant procedure works at 
the level of the action of the theory. It relies on the red- 
erivation of the field equations from an action defined 
within a closed volume of space time: when one follows 
the usual procedure and applies the Gauss-Stokes theo- 
rem, this confinement leads to finite terms (as one can 
no longer evoke that these vanish at infinity). In order 
to obtain the same equations of motion, these undesired 
quantities must be countermanded by a suitable bound- 
ary contribution, defined solely on the hypersurface — 
the Gibbons- York-Hawking term [36] (see Ref. [35] for a 
general derivation in the context of f(R) models). Fi- 
nally, by varying the action with respect to the induced 
metric on the later, one straightforwardly obtains the 
sought for second junction condition. 

In this study, one may quickly conclude as to what 
method is more profitable: if one adopted the first pro- 
cedure, the presence of the term A /il/ £ in Eq. (9) would 
lead to the appearance of derivatives of the Dirac delta, 
which are functionally defined as 6'(l)f(l) = —8f'(l). 
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This would lead to extremely cumbersome calculations, 
thus favoring the second approach outlined above. 

With the above in mind, one varies Eq. (1) with the 
adopted forms for fi = R and ji = 1 + eR/n and the 
constraint <5g M „ = in the hypersurface dV. Using 



SRa^v^srZp-g^srL] 



(36) 



one gets 



SS = J - s /~^gd i x x 



(37) 



k + -C^ R af 3 - ^nRg a /3 - ^ fl + e— ) 



H 1 



The second line in the integral above, dubbed SS2, may 
be integrated via the Gauss-Stokes theorem. For this, 

a/3' 



one first uses the definition of r' y g , obtaining 



U ) v i- k j v 

{-(V a V°£)g a pSg a P\ - { -(V a V p C)6g a ^} = 
Ik J y I k J y 

{n a ( K+ ic)g a ^(Sg^)} gv - 
{n a ( K+ ^)^(Sg^)} dv + 

+ {^(V a C)Sg^} dv - {^C)g aP Sg^} dv + 

\-(V a V°£)g af) Sg a P} - { -(V a V p C)5g^\ 
U J y I k > v 

{n a (k + ^C) g a ^ a (6g aP )} gv - 

{n„(. + i£)v,(^)} w + 

{-(V a V°C)g a p5g^X - {-{V a V p C)5g^\ 
U J y I k Ji 

having used 5g a P = at the boundary. 
Thus, the variation of the full action reads 



g a0 n(6g aP )-V a V p (5g af3 ) 



(38) 



For convenience, the integral symbols are substituted 
by an abbreviated notation, 



/ Xy/=g~d i x = {X} 
Jv 

I X\^hd 3 x = {X} 

JdV 



(39) 



where h is the determinant of the induced metric h a b] 
with this notation, the Gauss-Stokes theorem reads 



{V M X} y = {n^X} 



(40) 



Integrating by parts and twice using this theorem, one 
gets 



ss 2 = {(« + ^c) ^y ap sri ft -srsr^]} v = (4i) 

K+-£)g aP \7 a \7-(Sg a P)} - 
K 1 J V 

K+-£)V a Vp(5g a P)} = 
hi / J y 

K+-C)g a0 V°{5g a ( 3 )\) - 

K I i J V 

-(V„£)V^(^)) - f-(V CT £) ffQ/3 V ff (^)} = 
k J y I ft J y 

u CT («+-£) <7 Q /?V CT (<^)) 

V K J J i 



if R\ e . „ 
- 1 + e— T Q(3 A Q(3 £ 

Z \ ft I ft 



(42) 



{ (ft + lc) [n a9a ^(Sg^) - n a V (Sg^)] } 



The first two lines yield the modified field Eq. (9); the 
last line must be balanced by an adequate boundary 
term, given by 



6S b = 6{(n + -Cj [n a d (dg a ^ 



n<jg a f3 



d"{5g^)}} 



dV 

(43) 



since on the boundary Sg al3 = and thus V r7 {8g a ^ 3 ) 

dASg af3 )- 

Since Sg a ^ does not vary in the hypersurface, its 
derivative along the tangent vectors vanishes, 

<%4K*dfi(6g*a) = , (44) 
so that, from Eq. (34), one obtains the simplification 



( g <*0 -n a n? 3 )df > (6g aa ) = 
d a (6g aa ) = n a n?d p (6g aa ) 



(45) 



implying that 



dv 



n a dp(5g a P)-n a g a pd°(5g a P) 
n a (n a nfj - g a p)d° \5g af} ) = 
-n° ' e a a e\h ah d a ^g afi ) . 



(46) 
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As the variation of the trace of the extrinsic curvature 
K = h ab K ab with respect to g a/3 is given by 



1 



SK = -n° e a a e%h ab d a (6g afi ) 



one gets 



(47) 



one may now use Eq. (52) to fix the mass M, related 
to the Schwarzschild radius R s = 2GM; the former is 
expected to differ from the gravitational mass, defined as 
M = (4n/3)pRl 



VII. LAGRANGEAN DENSITY CHOICE £ = p 



5S b = -2{(k+ -£) 5K} 



dV 



(48) 



Since 28C = (g a pC — T a p)5g al3 =0on the boundary, one 
finally obtains the full action including boundary terms 
(written again in standard notation), 



S = 









I 


nR + 











<)Y 



K H — C ) K 

K 



\/—gd 4 x — 
eanx(h abl £) 



(49) 



-hd A x . 



The term x(h ab , C) (with dimensions of mass) stems from 
the fact that one can supplement the boundary terms 
above with additional contributions involving h ab and C 
only, as these do not show up when variation with respect 
to g a p is performed (again, since SC = on the bound- 
ary). It is factored by ea as one expects it to be present 
only when these quantities are non- vanishing (i.e. one 
has a non- minimal coupling with C ^ p = 0). 

The second junction condition may be obtained by 
variation of the above expression with respect to h ab on 
both sides of the boundary; considering that there is no 
surface energy-momentum tensor S ab describing a bound- 
ary layer, one has 



S ab = - = sc = l^ ab 5h ab . (50) 

so that, after manipulating the tensors, one obtains 



K b = (l - -^P) Kb + ^lPK-ha b + (51) 
+ea [X ab + h ab (x - X)] , 

defining X ab = Sx/Sh ab and its trace X = h ah X ab . 

Computing the extrinsic curvature tensor for the inner 
and outer metric, 



, 



(52) 



+0 _ 



GM + RlR* 



KJ" = K 



R*JR*[R* - 2GM + R*(R)l] 



R*[R* — 2GM + R*(R)l] 

m 



The choice of Lagrangcan density C = p — >= merely 
leads to a shift in the definition of the cycloid parameter 
T] Eq. (24) via the spatial curvature k, Eq. (19): thus, 
the gravitational collapse is dynamically equivalent to the 
case of GR, with the fundamental difference that is does 
not lead to a singularity of infinite density, as seen from 
Eq. (16). 

Solving for the original time t, one obtains 



dt 
drj 



a 

Vk 



COS 7] 



2Vk 



T) + sin r\ 
2y/k 



Gravitational collapse ends at a time rjf 
translates into the usual expression 



2y/k 



6 k 
po 



(53) 
7r, which 

(54) 



For a positive coupling strength e, one finds that the final 
state of finite density pj is attained earlier than in OS 
collapse, since the spatial curvature is larger. 



A. Apparent and event horizon 

Following the preceding section, if£ = p— >cn = one 
expects that the trapped surfaces, apparent and event 
horizon all occur analogously to OS collapse, with the 
non-minimal coupling manifesting itself merely through 
the shifted spatial curvature fc, Eq. (19). Indeed, by 

introducing the new radial coordinate \ = arcsin ( Vkr, 

one sees that the metric Eq. (11) becomes 



ds v 



■ drf 2 + dx 2 + sin 2 x(d0 2 + sin 2 



(55) 



so that radial photons follow null geodesies which are 
straight lines in the (rj, x) plane, as in GR. 

The calculation of trapped surfaces and the ensuing ap- 
parent and event horizon proceeds accordingly: in par- 
ticular, the apparent horizon crosses the surface of the 
star when it has collapsed below the Schwarzschild ra- 
dius, and becomes fixed at this value. 

In order to determine the latter, one resorts to the 
second junction condition, Eq. (52): setting a = 0, one 
sees that the continuity relation [K ab ] — is recovered. 
In particular, using Eq. (52), the tt component together 
with the identification i?* = ar* yields 
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= 



GM + RlR* 



(56) 



R^R*[R* - 2GM + R*{R)l\ 
R s = 2GM = -2a 2 arl = a[k + (a) 2 ]rj = krl 



having used Eqs. (17) and (22) with a = 0; the same 
result of course arises from [K$g] = 0. Inserting Eq. (19), 
one gets 



M = 



-in 



Po 



M 



3 1 - eo/2 * 1 - e /2 



(57) 



showing that the mass of the spherical body, as inferred 
by an outer observer, is increased due to the presence of 
the non-minimal coupling. 

Remarkably, this result shows that a non-minimal cou- 
pling can break the no hair theorem: indeed, two stars 
with the same gravitational mass but different sizes (i.e. 
initial densities po) will produce black holes with unequal 
event horizons. 



VIII. LAGRANGEAN DENSITY CHOICE C = -p 

The choice of Lagrangcan density C = p — >= leads to 
a more convoluted effect of the non-minimal coupling, as 
the coordinate time t is related to the cycloid parameter 
r\ through 



t 



1 /(l + cos 77) [(1 + cost;) 3 + 16e ] 



2v / A?V (l + cos?;) 2 +2(3 + cos77)e 



drj . (58) 



Since the above cannot be solved analytically, one re- 
sorts to a numerical integration, yielding the relation 
rj(t) depicted in Fig. 1. Substituting r)(t) into a(rj) = 
(1 + cos?7)/2 leads to the the modified evolution of the 
scale factor for different values of eoj shown in Fig. 2; 
notice that a large deviation of r](t) with respect to its 
GR counterpart arises even for small values of eo (where 
the shift of the spatial curvature Eq. (19) is negligible), 
due to the additional term in Eq. (58). 

The relative decrease of the collapse time t / compared 
to the elapsed period i/os = Tt/2\/k ~ for the OS scenario 
is depicted in Fig. 3, with the former being given by the 
equality 77 (tf) = ir: notice that for large eo, one gets 



1 



cos 77 



drj = 



(59) 



/o v ft v 3 + cos 77 ' V2k 
leading to the displayed asymptotic behaviour. 

A. Matching with outer solution 

In the present case C = —p, one unfortunately finds 
that the inner FRW metric Eq. (11) cannot be suitably 




FIG. 1: Evolution of the cycloid time for different values of 
eo = 10~ 3 (full), e = 1(T 2 (small dash), e = 1CT 1 (medium 
dash) and e = 1 (large dash); dotted indicates e = 0). 




FIG. 2: Evolution of the scale factor for different values of 
e = 1(T 3 (full), e = 1(T 2 (small dash), e = 1CT 1 (medium 
dash) and e = 1 (large dash); dotted indicates e = 0). 




FIG. 3: Relative decrease of the collapse time t e /t e os — 1 
for different values of eo; notice the asymptotic behaviour 
t e -> n/{2Vk). 
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embedded with the outer Schwarzschild metric Eq. (27). 
To ascertain this, one resorts to the junction conditions 
uncovered before, namely the coordinate matching at the 
boundary i?* = ar» and Eq. (52), here repeated with 
a = 1: 



K b = 



1_ -^P) K ab + -^P K h ab + 

■e [X ab + h ab (x - X)} , 



(60) 



By evaluating the above condition, one should be able 
to read the constant M. However, this turns out to be 
unattainable, both for x = as well as for a number of 
candidates for this extra term to the boundary action Eq. 
(49), e.g. i = pon = h ab \7 a \7 b p. 

Conversely, one may attempt to solve the above for x, 
thus obtaining the additional boundary term that must 
be considered so that the constant M is recovered: al- 
though this is in principle possible, inspection shows it 
to be extremely cumbersome, with the foreseeable result 
producing an extremely convoluted and unfounded ex- 
pression on p. 



1. Comparison with OS collapse 

In the absence of the non-minimal coupling, there is a 
much more straightforward way to approach the problem, 
which indeed produces the same results as the painstak- 
ing derivation of the Gibbons- York-Hawking boundary 
action [36] and the ensuing second junction condition 

[K ab ] = 0. 

In the standard OS collapse, the absence of pressure 
indicates that the dust particles on the surface of the 
spherical body are free-falling along radial geodesies of 
the outer Schwarzschild metric, so that 



R* = y(l + COS7/) , 



(61) 



R 3 

— l —(n + sin 7/) 
8M v/ ' ' 



where rj is a cycloid parameter related to the proper time 
t of the infalling observer [37]; the latter is identical with 
the coordinate time of the FRW metric, r = t. 

Recalling the solution Eq. (26) and the definition of 
the inner cycloid parameter Eq. (24) (here repeated for 
convenience) , 



a(rj) 



1 + cos(?y) 



(62) 



2aep 



y/k V ° 3 + a ( a + l) ae 



(e.g. the first junction condition) is only valid for all 
times in the a = case (see Eq. (53)), with Ri = r\ 
(since a(0) = 1) and 



V = V 



R 



8GM 2Vk 



(63) 



The later leads to the result of Eq. (56), 2GM = krl, 
valid both for GR (with M = Mo) as well as the non- 
minimally coupled a = case. 

Following this approach, one traces the impossibility 
of recovering an expression for M when a = 1 to the 
mismatch between the definitions of the cycloid param- 
eters Eq. (24) and Eq. (61). In its turn, this signals a 
fault in one of the assumptions of the procedure depicted 
above — namely that dust particles on the surface of the 
spherical body free-fall according to radial geodesies of 
the outer metric. 

Indeed, in GR this stems from the condition of van- 
ishing pressure, p = 0; in the non-minimally coupled 
scenario, an effective pressure arises, as the rr compo- 
nent of the modified field equations does not vanish, 
Peff = 2Kg rr G rr =/= 0, for a = 1 — and as a result 
dust particles in the surface experience an additional 
force that displaces them with respect to outer radial 
geodesies. It is null for a = 0, so that the above dis- 
cussion is valid, as attested by the matching between the 
inner and outer cycloid parameters. 

This is more than a simple mathematical curiosity of 
the scenario under scrutiny, as it recalls a similar prob- 
lem in GR: the impossibility of matching the inner and 
outer spacetimes in the case of a gravitational collapse 
of a homogeneous sphere p = pit) with non-vanishing 
pressure p 7^ 0. 

In GR, this can be alleviated by the inclusion of a 
suitable boundary layer, i.e. a finite surface energy- 
momentum tensor S ab , as given by Eq. (50). Such a 
procedure may also prove helpful in the present context, 
although it shall not be pursued here: as it stands, the 
inability to suitably enforce the required matching shows 
that the gravitational collapse of a linearly minimally 
coupled homogeneous sphere is only possible if the La- 
grangean density of a perfect fluid is given by C = p, not 
£ = -p. 



IX. DISCUSSION AND OUTLOOK 

In this work we have described the dynamics of grav- 
itational collapse of a dust sphere under the influence 
of a linear non-minimal coupling, thus extending the fa- 
miliar Oppenheimer-Snyder collapse. We have examined 
the different effects that arise due to the choice of La- 
grangean density of matter, namely the use of C = — p 
or C = p. We main results are threefold: 



one finds that the relation i?„ = ar t , stemming from the 
continuity of the induced metric h ab across the boundary 



- In the C = —p scenario, the dynamics of gravita- 
tional collapse deviates from GR, due to the more 
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evolved dynamics. However, the usual dependence 
of the density on the scale factor p ~ a~ 3 remains, 
and point-like singularity with infinite density are 
attained. 

The C = p case is much more interesting: although 
the dynamics are qualitatively the same as in GR 
(with all relevant quantities shifted but the same 
governing equations), the energy-momentum ten- 
sor is not conserved. This leads to a modified de- 
pendence for the density and, as a result, a geo- 
metric point-like singularity (i.e. where the scalar 
curvature diverges) is attained with a finite density! 
Depending on the value of e, this can fall within the 
classical domain, p — > pf <C Mp, thus avoiding the 
need for a description of the quantum regime. 

- Analogously to the well-known Gibbons- York- 
Hawking boundary terms, we have found that an 
additional contribution to the action functional on 
the hypersurface signaling the surface of the spher- 
ical body must be considered. Its Lagrangean den- 
sity is of the form Cgy = (l + eC/n 4 ) K, with pos- 
sible, undetermined, additional terms depending on 
the induced metric and the Lagrangean density of 
matter. 

By varying these boundary terms with respect to 
the former, we showed that the extrinsic curvature 
is in general discontinuous across the boundary of 
the spherical body. 

- The interior description of the gravitational col- 
lapse in the C = p case is suitably embedded into 
the surrounding Schwarzschild spacetime via the 
continuity of the induced metric and the extrinsic 
curvature. This leads to a shift of the mass M 
of the spherical body (given by the gravitational 
potential away from it) with respect to the gravita- 
tional mass M : this modification depends not only 
on the coupling strength e, but also on the value of 



the initial density po: as a result, different event 
horizons arise after collapse for stars with the same 
initial mass, but distinct radius — thus breaking 
the no-hair theorem. 

The scenario C = — p is not so well-behaved: the 
matching of the inner and outer spacetimes turns 
out to be unfeasible, unless highly unnatural, ap- 
parently arbitrary extra terms are added to the 
boundary action. From a physical point of view, 
this can be related to the non-vanishing effective 
pressure that arises due to the non-minimal cou- 
pling — thus recalling the similar matching prob- 
lem found in the gravitational collapse of a homo- 
geneous sphere with pressure in GR. 

Future work should clearly address the issue of consid- 
ering the collapse of a non-homogeneous sphere, as well 
as one endowed with initial angular momentum and/or 
charge. The consideration of a thin shell could also shed 
further insight as to the issues found with the spacetime 
matching in the £ = — p scenario, as well as providing a 
test for the validity of Birkhoff 's theorem. 

A more evolved, non-linear form for the non-minimal 
coupling will undoubtably also provide for a more evolved 
phenomenology. In particular, we speculate that such 
generalization could lead to a rather interesting scenario, 
where the spherical body does not collapse to a vanish- 
ing size, but instead asymptotically attains a final, finite 
radius. 
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